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Abstract 

We compute the QCD form factor resumming threshold logarithms in B ^ Xc + l + vi 
decays to next-to-leading logarithmic approximation. We present an interpolation formula 
including soft as well as collinear effects softened by the non- vanishing charm mass which, 
in the A^-space, reads: 
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ml) I (Q2 - m2), /5 = m2/(Q2 _ ^2) q = + \pxX The 
function A{p; as) describes soft radiation collinearly enhanced off the charm quark and 
reduces to the standard double-logarithmic function A{as) in the massless limit /O — > 0. 
The function B{as) describes hard collinear emission off the charm quark, being the 
standard jet function. D{as) and A(as) describe soft radiation not collinearly enhanced 
off the beauty and the charm quark respectively. 



1 Introduction and Summary of Results 



Distributions in semileptonic B decays 
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such as for example the hadron-mass event fraction, often receive logarithmic contributions of 
infrared origin of the form 
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where mx^ is the final hadron mass and as is the QCD coupling evaluated at the hard scale. 
Even in weak coupling regime as ^ 1 (pcrturbative phase), the terms above tend to spoil the 
convergence of the pcrturbative series in the threshold region, the latter being defined as the 
one having parametrically 

mx, -C mi,. (3) 

The large infrared logarithms arc a "remnant" of an incomplete cancellation between real 
corrections to the decay distributions, integrated on the border of the phase space, and virtual 
ones. An all-order resummation is needed in order to have a reliable description of the spectra 
in the whole kinematical range. The infrared logarithms coming from soft-gluon emission are 
of the form 
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and formally diverge for mx^ — > m'^. Since there is at most one soft logarithm for each gluon 
emission, i.e. for each power of as, it follows that k < n in eq. (2). The infrared logarithms of 
coUinear origin. 
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never become infinite but may become large in the threshold region (3). Since there are at 
most two infrared logarithms for each power oi as , k + I < 2 n in eq. (2). 

Our main result is the following expression for the resummed QCD form factor in A^-space 
(also called moment or Mellin space): 
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We have defined the hadron variable with unitary range 

2 2 

y = 7# — h 7 

Q'' - mi 

which equals zero in the Born kinematics, and the mass-correction parameter 

(8) 
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where Q, the hard scale in the heavy flavor decay, is given by 

Q = Ex^-r IpxJ , (9) 

with Ex^ and being the total energy and 3-momentum of the final hadron state X^. The 
function y4(p; as) has an expansion in power of q;^, 

oo 

A{p-as) = (10) 

n=\ 

and describes soft radiation collincarly enhanced off the charm quark; it reduces to the standard 
double- logarithmic function A{as) in the massless limit: 

A{p-as) ^ A{as) for p ^ 0+ . (11) 

The first-order coefficient is: 

AW(p) ^ ^(l + 2p), (12) 
vr 

where Cp = {N^ - l)/{2Nc) = 4/3 for A^^ = 3. The second-order coefficient A^'^\p) (as weU 
as the third-order one A'^^^p)) is only known in the massless limit (see later). The function 
B{as) is the standard massless jet-function, describing hard collincar emission off the charm 
quark, while D{as) and A{as) describe soft radiation not collincarly enhanced off the beauty 
and charm quarks respectively. All these functions have a perturbative expansion in as^: 

oo oo oo 

Bias) = Yl ^^"^ «S ; D (as) - J] D^^^ a^s ; A {as) = ^^"^ «5 ■ (13) 

n=l n=l n=l 

The 0{as) coefficients read: 

^(i)^_3CV. A« = -^. (14) 

4 TT TT TT 

The form factor (6) aims at describing the three different dynamical regions in the semileptonic 
decay (1), which are identified by the value of the (ordinary) charm velocity Uc = Pc/Ec — 
without any generality loss, we can work in the beauty rest frame. It holds: 



These regions arc: 



^For a compilation of the massless coefHcients, see for example [1], while for the heavy flavor case [2]. 
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1. very slow charm quark: 



Uc > or, equivalently, p ^ 1 . 



(16) 



For Uc = there is no soft-gluon emission to any order in perturbation theory [3] ^. That 
is a coherence effect, namely destructive interference between emission off the initial and 
the final state. The physical picture is the following. The b decays into the c at rest 
and with the same color state of the b. Soft gluons "view" the b as well as the c as 
static color charges, and therefore do not see any acceleration or color-spin fiip, i.e. any 
change occurring in the decay. As a consequence, there is no radiation. Furthermore, for 
small Uc, radiation is proportional to u^, i.e. first order terms 0{uc) vanish [4]. Radiative 
corrections to our form factor vanish in the no-recoil point: 
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implying that this region is correctly described by the form factor; 
2. non-relativistic charm quark. Schematically: 
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The final state Xc contains the final charm quark together with soft gluons emitted at 
any angle with respect to it, because there is a significant soft emission. Since the QCD 
matrix elements do not have any collinear enhancement, Xc does not have a jet structure. 
In the threshold region y <^ 1, so that y <^ p and the form factor (6) simplifies into: 



crs,N{p,Q'^) = exp 
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3. fast charm quark: 

Uc<l or p < 1 . (20) 

Soft gluons are mostly radiated at small angle with respect to the charm quark and there 
is a jet structure of the final state. For p — > 0+ the form factor (6) reduces to the standard 
massless expression [2,5-7]: 
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"'''A[as{k')]^ + D[as{Q'y')] + 



By "soft" we mean a gluon with energy Eg < mc,mi,. 
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+ B[as{Q'y)] 



(21) 



where: 

Q ^ Ex + \px\ ^ '^Ex, (22) 



Ex - \px\ ^ rnjc 
Ex+\px\ ~ ^El 



where on the last members we have expanded for mx <^ Ex- The inclusion of first- 
order corrections in p on the r.h.s. of eq. (21) basically amounts to take into account the 
dead-cone effect in gluon radiation off the charm [8]. The form factor factorizes into the 
product of the form factor for a massless charm times the universal mass-jet correction: 

a^^(p, Q^) ^ ^T^^(0, Q^) S^ip, Q^) for p < 1 , (23) 

where ^ 

S^ip, Q') = exp / dy^-^ -{ [a {kl)] -B[a {pQ'y)] + 

Jo y [ JpQ^y'^ 

+ D[a(pgV)]|- (24) 

Let us note that in [8] the universal correction 6isf{p, Q^) was found for p <^ 1, while in 
this work we consider a specific process for any p. 



Since the ratio mc/m^ ^ 1/3^ 1/4 in the real world, i.e. it is not so small, region 3. is tiny. 
With typical values of the on-shell masses, mj, = 4.7 GcV and rric = rrib — uib + mo ~ 1.29 
GeV [9], one has a maximal charm velocity Mcmax — 0.86 corresponding to a Lorentz factor 
7c max — 2, or 0.081 < p < oo. With hadron kinematics, i.e. with = uib and rric = mr), one 
obtains instead: Wcmax = 0.78, 7c max — 1-6 and pmin — 0.14. 

Let us comment further upon the properties of the form factor. 

Since the y and A;^ integrals on the r.h.s. of eq. (6) extend down to zero, one hits the infrared 
singularity in the QCD couplings as{k'^) and as{- ■ ■y'^) — the well-known Landau pole: some 
prescription is needed in order to render the integrands well behaved. As it is usually the case 
with resummation formulae, eq. (6), as it stands, has an "algebraic" sense: upon expansion in 
powers of as, it allows to predict the logarithmic corrections which would be explicitly found 
in higher-order Feynman diagram computations. 

^The variable p used in this work and the variable r = m-^/Q^ = + p) used in [8] coincide to first order. 
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In general, we expect less radiation to be emitted in the decay (1), because of the rather 
large charm mass, compared to the charmless channel 



B ^ X^ + l + ui. (25) 

As a consequence, the typical Sudakov effects, namely suppression of non-radiative channels 
and broadening of sharp structures, are expected to be less pronounced for our process. In 
principle, for the decay (1), one has single-logarithmic corrections, which arc not strong enough 
to shift the peak of tree-level distributions. For y ^ they produce indeed a form factor of 
the form 

a(y) ^ -^e-^^'^^-y ^ (26) 

which still has an (infinite) peak in 7/ = 0, like the lowest-order form factor a^^\y) — 5{y). 
That is to be contrasted with a double-logarithmic form factor ^ d/dyexp{—as log^y). 

It is remarkable that the limits of the A;^-integral on the r.h.s. of eq. (6), Q^y^ jiX + p) and 
Q^y^ I {y+p)-, are equal to the arguments of the QCD couplings entering the "sub-leading" terms 
(i.e. those without the explicit k'^ integration); that seems a general feature of resummation 
formulae. A consequence of this fact is that only 3 out of the 4 resummation functions A, B, 
D and A entering ctn are actually independent (see sec. 3.3). Similar considerations can be 
repeated for the soft form factor in eq. (19), in which the scales Q^y^ j (1+p) and Q^y'^ j p appear 
as limits of the integral as well as arguments of as in the D and A terms. In the latter case, 
only 2 out of the 3 functions A, D and A are independent. 

By looking at the ^/-integration on the r.h.s. of eq. (6), one finds that the threshold region 
cannot be described within a "pure" perturbative framework when the QCD coupling comes 
close to the infrared (Landau) singularity, namely when 



KcD , (27) 



1 + p -^'2^^' 

where Aqcd is the QCD scale, i.e. when the final hadron mass becomes as small as ^ 

+ AqcD^Q^-ml. (28) 



m 



~ m" 

NP 



These effects are related to soft interactions only, because hard collinear terms (oc B{as)) 
are controlled by a smaller coupling, evaluated at the larger scale Q'^y'^ /{y + p). These soft 
interactions are non-perturbativc in region (28), and can be factorized into the shape function 
for a massive final quark [10]. Note that for rric — > 0, the slice (28) reduces to the well-known 



"^For Q = m,}, and Aqcd = 300 MeV one obtains for the current values of the b and c masses mxjjyp — 1.7 
GeV. 
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one ~ ^qcdQ ■ In region (28), also soft interactions between the b quark and the hght 
degrees of freedom in the 5-meson, namely the valence quark, are important. That is the 
well-known "Fermi-motion" , characterized by momentum exchanges A;^ in the i?-meson of the 
order of the hadronic scale, 

I/CmI ~ ^QCD. (29) 

Physical intuition would suggest that Fermi motion is related to the initial S-meson state only, 
i.e. that it is independent on the mass of the final quark. In quantum field theory that is no 
true and it is not easy to find a general relation between cq. (21) with the B term dropped, 
describing soft interactions for rric — > 0, and eq. (19), describing soft interactions for rric ^ 0. 

It is technically simpler to compute first the form factor for the radiative case 

S ^ X, + 7 (30) 

for rris 7^ and then to extend the result to the semileptonic case (1), rather than to deal 
directly with the latter process. That is the strategy we will follow. 

The plan of this note is the following. 

In sec. 2 we consider soft effects only, which are in some sense "leading", as they become 
formally infinite at the border of the phase space, for the radiative decay (30). 

In sec. 3 we include hard coUinear effects, which never become infinite but are logarithmically 
enhanced for small strange mass, as discussed at the beginning of this section. 

In sec. 4 we derive the resummed form factor for the process of main physical interest, namely 
the semileptonic decay (1), by means of a generalization of the radiative case. 

Finally, in sec. 5 we draw our conclusions and we give an outlook to possible developments of 
the results, both of theoretical and phenomenological nature. 



2 Soft-Gluon Resummation 

In this section we consider the resummation of threshold logarithms coming from multiple soft- 
gluon emission to all orders of perturbation theory. The QCD form factor factorizing such 
terms has a perturbative expansion of the form: 

00 n—l k \ 

n=i k=o \ y J + 
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S{y) + ^""^ 



TT 



(l + 2p)log(i±^) -2j +0{al). (31) 



As discussed in the introduction, we first consider the simpler radiative decay (30), for which 
y is naturally defined as 



2 2 



y ^ (32) 



and the mass-correction parameter p as 



ml -ml 



ml 



P = 2 ' 2 ■ (33) 
— mg 

The c„fe(p)'s are functions of p and the plus distributions are defined as usual as: 

P{y)+= lim \e{y-e)P{y)- S{y - e) f dy'P{y') ] . (34) 

Let us note that we do not have here a double-logarithmic problem as in the massless case 
(p = 0), but a single-logarithmic one: there is at most one infrared logarithm for each power 
of as- The additional singularities for p — > 0"*" are "hidden" in the coefficients c„fc(p)'s, which 
indeed diverge logarithmically in that limit — see the last member in eq. (31). Note that 
the square bracket on the last member of eq. (31) vanishes for p ^ oo (see region 1. in the 
introduction). By leading logarithmic soft (LLs) approximation, we mean the resummation of 
all the terms in eq. (31) with A; = n — 1, by next-to-leading logarithmic {NLLs) approximation 
also those with k — n — 2 and so on. 

2.1 Single-Gluon Emission 

In this section we consider the emission of a real gluon, i.e. the process: 

b ^ s + g + -f. (35) 
We will recover the soft contributions factorized in as above. 

2.1.1 Kinematics 

The energy of the strange quark in the tree-level process 

6 ^ s + 7 (36) 
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is fixed to 



and its 3-momentum to 



= (37) 
2 1+p ^ ' 

^(0) = 1^(0)1 _ e;(0) _ .goN 
Vs - Ws I - ^7 - 2(1 ^^^^ 

The strange quark and the photon are emitted back to back in the decay of the h quark at rest. 
Note that Ef^ + = . Let us now consider kinematics of single-gluon emission. The 
maximum value of the gluon energy is 

^2(1+/)) ^ ^ 

and it corresponds to the configuration with the strange quark and the gluon back to back with 
an accompanying soft photon. Let us define the unitary gluon energy lo as: 

IE 

u;=^(l + p) (40) 

and the unitary angular variable t as 

1 — cos Q / 1 ^ 2 



t = ^ j for e « 1, (41) 

where 9 is the emission angle of the gluon with respect to the final (on-shell) strange quark. 
Let us now restrict our attention to kinematics in the soft limit. Eg and Ps are only mildly 
modified by the emission of a soft gluon, so that we can neglect recoil effects and set 

Es^Ef^; Ps^pf. (42) 

Furthermore, the direction of motion of the strange quark is modified by gluon emission in 
negligible way. In terms of the variables defined above, we can write 

y = (43) 

^ 1+p ^ ' 

Note the linear relation between y and a;, which is a consequence of the soft limit, and that 

y ^ <^ u; ^ 0. (44) 

That implies that the threshold region y <^ 1 selects kinematical configurations which are 
dynamically enhanced by the soft divergencies. At fixed one has for the angular variable t 
the whole range < t < 1. Finally, note that for any hadron final state: 

Ex, + \px,\ = mi,. (45) 
That is a consequence of the massless "probe", namely the photon (g^ = with q = p^). 
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2.1.2 Eikonal Current 



In the soft limit a; <C 1, the emission of a real gluon off the b and the s quarks is described by 
the eikonal current: 

where Tb and Tg are the color generators of the b and s quarks respectively [8] . By squaring 
the eikonal current, integrating the phase space over the azimuth and summing/averaging over 
helicities and polarizations, we have the following contribution to the differential form factor 
from the emission of a real gluon: 

Let us briefly comment upon the result above. The soft singularity for a; — >^ of the QCD 
matrix element squared is screened by the kinematical constraint, which is infrared safe {y ^ 
=^ a; 7^ 0) [11]. Also the coUinear singularity is screened by the kinematical constraint, 
which is collinear safe (?/ 7^ ^ t 7^ for p = 0), as well as by the strange mass, because 
the expression (47) diverges for t —>■ —p only, which is outside the physical domain {t > 0) for 
p > 0. 



2.2 Inclusive Gluon-Decay Effects 

In the previous section we considered the 0{as) correction to the tree-level decay (36), namely 
the emission of a single gluon. Higher order corrections in as basically involve multiple emissions 
of the following two kinds: 

1. primary emissions, i.e. direct emissions of gluons directly attached to the colored 6 — > s 
line; 

2. secondary emissions, produced by the splitting of a gluon into a gg, qq pair, or by gluon 
emission off secondary quarks or antiquarks, etc. 

As we are going to explicitly show in sec. (2.4), the main effect of primary emissions is to 
exponentiate the single-gluon distribution. Let us consider in this section secondary emissions, 
namely the "decay" of a primary (virtual) gluon with momentum pg into secondary partons. 
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The amplitude for the emission of one gluon off the strange quark is enhanced when the 
scalar part of the s-quark propagator, 

= (48) 

becomes large, i.e. when the virtual s quark is close to its mass-shell. This denominator reads 
in the case of the emission of one virtual gluon with momentum pg and invariant mass squared 
ml^pl: 

^s{Ps+P,) = ! , ^ , .. T^ (49) 



where we have taken the strange quark after gluon emission to be on-shell, i.e. Pg — m^, since 
we are considering a single primary emission. On the last member we have made the soft 
approximations for Eg <^ Eg discussed above and we have dropped the term, which comes 
out to be negligible a posteriori (see later) . Since infrared singularities originate from quasi-real 
partons, i.e. from partons with a small virtuality, we can assume that <^ E'^, to obtain [12]: 

^s{Ps + p.) - — , (50) 



(l+p)2 mg 



where we have made the small-angle approximation 1 — 2f ~ 1 in the term multiplying m^, since 
the latter is already a correction. Now it comes the main point: the 0{as) correction involves 
a real gluon, having = 0, and is logarithmically enhanced in the infrared regions discussed 
in the previous section. In order not to spoil the lowest-order enhancement, the denominator 
has to remain small, i.e. the correction term in the square bracket has to be much smaller than 
one: 

where we have replaced the strong inequality with an ordinary inequality since we work within 
logarithmic accuracy. The condition above specifies the upper limit on the gluon mass which 
does not spoil the leading logarithmic structure of the cascade. It is immediate at this point to 
justify neglecting the term on the last member of eq. (49) : 

2 

m 

^ < a; < 1 . (52) 



2ps ■ Pg 

The amplitude for the emission of a soft gluon off the initial h quark (at rest) contains the 
scalar part of the b propagator 

^b{Pb - Pg) = 7- -T^ — 2 = 9^ W _L ^2 ' (^"^) 
\Pb-Pg) -^b -ZrUbl^yg + m^ 
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which has a different enhancement with respect to the s one; in particular, it involves only 
the gluon energy and not the angle with respect to the s quark. However, for ^ 0, the 
kinematical constraint of fixed y is no more linear in uu (cfr. eq. (43)) and reads: 



1 + P 



t + p + 



[1 + pfm: 



(54) 



In particular, for y — > O"*" one does not select anymore the limit a; — > 0"^ unless the gluon mass 
satisfies the constraint (51). As a consequence, for kinematical reasons, also amplitudes for soft 
emission off the b are enhanced only if condition (51) holds. 

The conclusion is that, as long as the restriction (51) is satisfied, one can sum over all 
possible gluon decay channels. The y-distribution is an example of "inclusive gluon decay" 
quantity. By including real as well as virtual gluon-splitting processes, the sum of cuts of the 
gluon propagator can be written as its discontinuity, so that the tree-level QCD coupling is 
replaced by the following momentum-dependent effective coupling [13]: 



as ^ as (A;') = — 



fc2 



dnig Disc 



1 



'"?''«"'8(Aa)j 



(55) 



where (3o = (11/3C^ — 2/3n/)/(47r) with Ca = = 3 and we have defined the transverse 



momentum squared as the largest possible m^: 



2 1 2 ^ ~^ P 

max rUg — m^uj 



(56) 

terms in gluon polarization 



(1+p)' 

On the last member of cq. (55), absorptive effects (the "— ivr 
function) have been neglected [14] ^ . Eq. (55) implies that, in order to include leading gluon- 
decay effects, one has to evaluate the running coupling at the transverse momentum squared. 
The latter has the range: 



ml y2 



^eK^^. (57) 
1 + P P 

Let us note that this range shrinks to zero for p — > oo, killing any coUinear enhancement. 
Finally, let us observe that in the massless limit for the strange quark, p = 0, one recovers the 
definition of transverse momentum given in [6] . 



2.3 Effective Single- Gluon Distribution 

In this section we include gluon branching effects, i.e. secondary emission, in the single-gluon 
distribution by means of the replacement (55). The real contribution to the form factor therefore 

^The over-all sign on the r.h.s. of eq. (11) of [14] should be minus. 
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reads: 



TT Jn Jo 



UJ 





r 2 2 1 
mi UJ T 






< 





UJT 



1 + P 



where we have changed variable from t to 

T = t + p. 

Let us consider the integration of the above terms in turn: 



(58) 
(59) 



1. after integrating over r and changing variable from a; to A;^, the first term in the square 
bracket on the r.h.s. of eq. (58) reads: 



1 Cf 



2 2 

"'by 



-^(l + 2p) / ^as{k'). 



(60) 



Note that for p » 1 the coupling becomes large and we leave the perturbative phase of 
QCD even outside the threshold region, i.e. for large y ~ 0(1); 

2. the second term on the square bracket (related to soft-gluon emission off the 6), after 
integrating over uj the 5-function, reads: 



1 Cp I , 
/ dr as 

y jp 



r 



( ml 
y TT \1 + p 



(61) 



3. the last term in the square bracket (related to soft gluon emission off the s quark) is 
treated in a similar manner and reads: 



Cfp{1 + p) f'+^dr fmy\ ^ 



TT 



y 



y TT 



as 



2 2 

mf y^ 



(62) 



Let us note that in 1. we have rewritten the logarithmic integration over r in terms of k^, while 
in 2. and in 3. we have explicitly computed the non-logarithmic integrals over r. 



By summing all three terms, we obtain: 

(^s^\y; p,ml) 



2 2 

1 1 dk' 



k^ 



' '''' A^^\p)as{k^) + Di^^as l^] + A«a. 



+ P 



(63) 

In the frozen coupling limit as (• • •) —>■ as or, equivalently, in first order, the above expression 
reduces to the real contribution in the last member of eq. (31): 



(^s^\y,P, as) 



CpOis 1 
TT y 



(1 + 2p) log 



1 + P 



- 2 



+ O {al) 



(64) 
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2.4 A/^-space and Exponentiation 



The exponentiation of threshold form factors in gauge theories is a consequence both of dy- 
namical and kinematical factorization properties [11,15,16]. In QED, multiple soft-photon 
matrix elements factorize into the product of single-photon matrix elements because of inde- 
pendent emission (Poisson distribution). In QCD that is no longer true: soft gluons are not 
emitted independently because of gluon splitting, but gluon correlations largely cancel in the 
form factor, leading to an "effective" matrix element factorization. On the kinematical side, to 
have exponentiation, the constraint for multiple soft emission must factorize into the product 
of single-gluon (or single-photon) kinematical constraints. The latter property does not hold 
in the physical, energy-momentum, space and, to have a consistent resummation, one has to 
transform the form factor to A'"-space: 



f\l-y)^-'as{y)dy. 
Jo 



(65) 



The QCD form factor in the original (physical) space is recovered by means of a standard 
inverse transform: 

1 /-C+ioo 

^s{y)^,^. dN{l-y)-''as,N, (66) 

where C is a real constant chosen so that all the singularities of as,N he to the left of the 
integration contour. 

Up to now we computed only real emission contributions to the form factor. Virtual correc- 
tions do not need to be explicitly evaluated: they can be included by imposing the vanishing 
of the first moment {N = 1) of the form factor, implying the replacement 



N-1 



N-1 



1. 



(67) 



By inserting the expression for o"^^^(y) obtained in the previous section and exponentiating 
the result in order to take into account multiple primary emissions (real and virtual ones), we 
obtain: 



crs,N{p,ml) = exp 



r'dy 

Jo y 



7V-1 




as 



2 2 ^ 



(68) 



Eq. (68) consistently resums soft logarithms in leading logarithmic (LLs) approximation ^, i.e. 



^The form of the resummation coefficients given in eq. (68) is not unique. An alternative form involves a 
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all the soft terms having the form in A'^-space 

a^log^'N. (70) 

For a complete next-to- leading logarithmic (NLLs) approximation, i.e. the inclusion of all the 
terms 

«glog"-^7V, (71) 
one has to include second-order corrections to the above first-order terms, i.e. 

A^^\p) as ^ A^'\p) as+A^^\p) 4 ; D^'^as ^ D^'^ag+D^^^al ; A^'^ag ^ A«q;5+A(2)q;| . 

(72) 

These higher-order corrections describe residual multi-parton interactions which are not taken 
into account neither by the exponentiation of the one-gluon distribution nor by the insertion of 
the effective gluon-jet coupling osik"^) ^. The second-order corrections are, as far as we know, 
not available at present in QCD ^. However, an improvement of eq. (68) can be obtained by 
using the two-loop approximation for the QCD coupling as and by adding the second-order 
correction A^'^^(p) in the limit p ^ 0, which is known since a long time [17-19]: 



7r2 



^ 36 2 7 18 ^ 



(73) 



where z{2) = tt^/G and n/ is the number of active quark flavors. Let us note that the term 
proportional to A^^^ gives a contribution on the r.h.s. of eq. (68) 



^(=>„|log(i±^), 



(74) 



which vanishes in the limit p — > oo, and therefore does not spoil the vanishing of the exponent 
of the form factor in that limit (see eq. (17)). Soft logarithms can formally be resummed to 
any accuracy by including higher order corrections to the above functions, namely: 

2 2 

{^b y 



coefficient A^^^ equal to the massless one and coefficients D^^^ and A^^^ depending on p: 

id) = £l, dW^p) = _^ (i_^logi±^) ; aW(p) = fi_plogi±^) . (69) 

TT TV \ p J TV \ p J 

Note that &^\p) and A(i)(p) vanish for p oo. 

^That is similar to the (more familiar) inclusion of higher-order corrections to a one-loop anomalous dimen- 
sion: 7^"^ as ^ 7 (as) = 7^°^ as + 7^"^^ q;| + ■ • ■ . 

^Virtual corrections to heavy flavor decay have been computed in two-loop order in soft approximation 
within an effective theory framework for n/ = (purely gluonic case) in [4]. The relation to QCD distribu- 
tions is straightforward to 0{as), but beyond leading order we do not fully understand the regularization- 
renormalization scheme effects. 
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The functions D {as) and A {as) coincide up to first order and, on the basis of physical argu- 
ments, it was conjectured in [8] that this property extends to higher orders, i.e. that ^ 

A {as) = D {as) . (76) 

We do not deal with this problem here because A^'^\p) is not available at present and therefore 
our analysis is limited to including D^^^ and A^^) only. 

The soft resummation formula in eq. (75) is the main result of this section. Let us com- 
ment upon it. The term proportional to A{p; as) contains a logarithmic integration over the 
gluon transverse momentum squared and becomes large for p small: it contains a kind of 
"p-evolution" , while in the massless case one has an evolution in y (cfr. eq. (86)). The contri- 
bution in D{as) comes from soft emission off the beauty quark, while that one in A.{as) from 
the strange quark. 



3 Resumming Final-Mass Effects 

In this section we derive a general resummation formula by using the dipole factorization 
formulae for the massive case derived in [20,21] (see also [22]). By taking the photon in the 
final state as the spectator parton, one obtains from eq. (5.16) of [21]: 



iv) = f 
y TT 

y+p 



l + 2p 1 + z p 



dz + 



-{l-y)z 2 y 
+ (soft emission contributions off the b, oc D^^^), (77) 

where z is the normalized projection of the strange 4-momentum along the (reference) photon 
4-momentum: 



{Ps+Pg)-Pj' ~ {Ps + Pg) ■ 



z,^ 'J' ''^ = 1-z,. (78) 



The factor 1 + 2p is not present in eq. (5.16) of [21]: we have inserted it by hand in order to 
take into account the specific soft structure of heavy fiavor decay, namely one heavy leg in the 



^An explicit check of the above relation requires a massive two-loop computation, which is outside the scope 
of the present paper. 
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initial state and one lighter leg in the final state; this process was not explicitly considered in 
[21]. ^° Let us now sketch the integrations of the above terms in turn: 

1. by keeping logarithmic terms only, the first term in the square bracket on the r.h.s. of 
eq. (77) reads: 

Cf«s 1 + 2p du ^ _ Cpas 1 + 2p / y + p \ ^^g. 

TT y Jy ^ TT y Vi + P/' 

where the unitary gluon energy reads in terms of light-cone variables: 

uj = l-{l-y)z^l-z + y. (80) 
Eq. (80) is easily proved by using the following expression for y: 

y = = ^111^ ; (81) 

Ps-P^ + Pg-P^ + Ps- Pg 

2. The second term on the r.h.s., related to hard coUinear emission off the s quark, reads: 

Cpas 



TT 



3 1 p 

+ T 



(82) 



4l/ + P ■ 4 (y + p)2 

The second term in the above equation is always smaller than the first one and it does not 
produce any logarithmic contributions. Therefore, it can safely be neglected. For y ^ p. 
we re-obtain the well-known massless contribution, — 3Cira5/(47ri/), while in the opposite 
regime, y <^ p, the above contribution can be neglected because there is no more a small 
denominator. Let us stress that all that is in complete agreement with physical intuition; 

3. The last term on the r.h.s., related to soft emission not coUinearly enhanced off the strange 
quark, reads: 



^ \y y + p 



3.1 Inclusive Gluon-Decay Effects 

The gluon transverse momentum squared in conveniently written as 



e^^, (84) 

l + p ^ ' 



us stress that 1 + 2p ^ 1 in the quasi-colhnear Hmit {k\_ — > 0, m ^ with k^jm const). That 
imphes that this factor cannot be determined with quasi- coUinear factorization but only with soft factorization 
because it is process dependent. 
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where to is expressed in terms of light-cone variables in eq. (80) (see also [23]). As in the soft 
case, we include gluon-branching effects by replacing the tree-level coupling with the effective 
coupling (55) in the one-gluon distribution (77): 



1. The double- logarithmic term on the l.h.s. of eq. (79) then reads: 

2p) 



1^(1 

y TT 



2 2 



2 2 T2"«5(K ) • 



(85) 



In the soft hmit (y — > at fixed p G (0, cxd)), the above expression reduces to the 
corresponding eikonal contribution derived in eq. (60) of the previous section: the only 
difference lies in the upper limit. In the massless limit (p — > with y ^ 0), the above 
formula reduces instead to the standard (massless) contribution: 



y Jmly^ 



(86) 



2. The determination of the argument of the QCD couphng entering the subleading terms 
goes as follows: 



1 Cf 



y 

y+p 

ICp (mly^ 

as\ — ■ — 

y TT \y + P 



2 y 



as 



mly{l- {l-y)z) 



1 + P 



dz 



p 

y+p 



1 + Z 



P]dz + 0l^ 



(87) 



where on the last member we have replaced the following expansion of the QCD coupling: 



as 



as 



1 1 

y + p 



y + p 



(i + p)y 



(88) 

The 2;-integral on the r.h.s. of eq. (87) is equal to the one in lowest order we have already 
calculated and the higher-order terms in the expansion (88) produce terms beyond our 
accuracy, namely OipCgjy). Let us remark that the rule (55) can be applied also to 
describe the branching of a hard gluon emitted coUinearly to the s-quark, because the 
enhancement is given again by the denominator of the s-quark propagator (see sec. 2.2). 



By explicitly adding the soft emission contribution off the 6 obtained in the previous section, 
the complete distribution reads: 



2 2 



2 2 



+ P 



1 

y 



(89) 
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as 



1 2 



as 



1 



y + pj y + p 



3.2 A/^- Space and Exponentiation 



By repeating the computation made in the soft case, one obtains: 



2 2 

a^ip; ml) = exp / dy[{l - yf-' If J A^'^ (p) as{k 



dk 



1.2 



+ 



y + p 



Ad) 



as 



2 2 

y + p 



as 



2 2 

mjy^ 
1 + P 



+ 



y + p 



5(1) 



as 



2 2 - 

ml r 

y + p, 



+ 



(90) 



What about the inclusion of higher orders? The first point is that there is a different counting 
in the soft region p> 0(1), which is a single logarithmic region, and in the collinear (p = 0) 
or quasi-collinear regions (p ^ 1), which are double logarithmic ones. For instance, for large 
p ~ 0(1), the term proportional to A^^\p) is of the same order of the Z^*^^) and A*^^) terms. On 
the contrary, for p ^ 1 the A^) term is leading because it is multiplied by a double log, while 
L)(i) and A^^) are subleading because they are multiplied by a single log. That complication 
arises because our formula is a smooth interpolation between two different dynamical regions. 
An improvement can be realized as in the soft case by means of the replacement: 

(p) as ^ A^') (p) as + A^^^ (0) 4 , (91) 

with A^'^\Q) given in eq. (73). The above substitution allows a complete NLL resummation in 
the massless case {p/y — > 0) and a partial NLL resummation in the soft case {y/ p — > 0). 

Infrared logarithms of subleading orders can be formally resummed by including higher- 
order corrections to the functions in the exponent, similarly to the soft case: 



(^N{p]m^) = exp / dy 

10 



y) 



N-l 



+ 



y + p 



1 



2 2 

y + p 



as 



2 2 " 

m^^y 



1 + 



as 



2 2 

y + p 



+ 



y + p 



B 



as 



2 2 

y + p 



+ 



.(92) 



3.3 A Symmetry of the Resummation Functions 



Eq. (92) is invariant under the simultaneous substitutions: 



A{p-as) A{p,as) + P{as)-^F{as) ] D{as) ^ D{as) + F{as); 
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A{as) ^ A{as) - F{as) ; 



B{as) ^ B{as) - F{as) , (93) 



where F{as) is an arbitrary function having a power-series expansion in the QCD couphng, 

oo 

F{as) = (94) 

n=l 

We have defined the QCD /3-function coefiicients with an over-all minus sign: 

= /?(«5) = -Poal - Pial (95) 

The above relation implies that only 3 out of the 4 functions entering eq. (92) are independent; in 
particular, one can choose F{as) in such a way as to make one out of the functions B{as) , D{as) 
and A.{as) vanishing, or to reduce the function A[p^as) to its lowest-order term A^^\p)as- 
The above symmetry is a generalization of that one of the massless formula (p = 0) involving 
only the functions A, B and D found in [24]. In NLL, the above substitutions read: 

(96) 

Finally, the rcsummation coefficients are also modified by a change of the renormalization 
scheme of = '(1 + Cas ' + •••)> this problem has been investigated in [25]. 



4 Semileptonic Decay 

This is the central section of the paper, in which we extend threshold resummation from the 
radiative decay to the semileptonic one. By means of general field-theoretic arguments, one 
can show that the hard scale Q in heavy fiavor decays is proportional to the total final hadron 
energy [2,7]: 

Q^2Ex, + 0{-^y (97) 

where f = s,c. The main difference between the radiative decay (g^ = 0, q = p^) and the 
semileptonic one (g^ > 0, q = Pi + Pv) is that in the latter case the lepton-neutrino pair can 
take away from the QCD final state a large fraction of the initial energy provided by the 6-quark 
mass. We may fix the power-suppressed terms in eq. (97) by imposing that the hard scale Q 
exactly coincides with the beauty mass when = (see eq. (45)), so that: 

Q ^ Exf + \pxf\ rrib when = . (98) 
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The hadronic variable y and the mass-correction parameter p must therefore be defined by 
replacing m;, with the general (g^-dependent) hard scale Q: 



'^'xs ~ "^s "^^c ~ "^c "^s "^c 

y 2 9" ^ 9" ) P 2 9 ^ /09 9 ■ 



One then obtains the resummed QCD form factor in A^-space given in eq. (6) of the introduction, 
the main result of this work. This equation systematically resums all threshold logarithms 
occurring in the semileptonic B decays for any value of the ratio mc/m,}, G [0, 1]. It interpolates 
in a smooth way dynamical regions 1., 2. and 3. described in the introduction. 

If the charm quark is not (too) fast, as it happens in most of the phase space of semi- 
leptonic b ^ c transitions, one can use the simplified resummation formula (19), which resums 
soft logarithms only. The latter correctly describes regions 1. and 2. but, unlike the full result 
(6), cannot account for the massless limit p —>■ 00, i.e. for region 3. In the soft limit: 

EiJ" = ; |p- | ^ |p (o.| = V^if'-W) , (joo) 

where A(a, b, c) = a"^ + b'^ + — 2ab — 2ac — 26c. It is straightforward to check that the above 
formulae reduce to the one for the radiative case in the limit — > 0. 



5 Conclusions 

We have presented in eq. (6) an expression for the QCD form factor resumming threshold 
logarithms in the semileptonic decay 

B ^ Xc + l + (101) 

The resummation is valid for any chosen value of the ratio iric/mi) G [0, 1] and reduces to the 
standard next-to-leading logarithmic approximation for itlc 0. Coherence effects, occurring 
for a small relative velocity of the c and b quarks and leading to destructive soft-gluon inter- 
ference, are correctly incorporated in the form factor, whose corrections vanish in the no-recoil 
point. Our form factor is basically a smooth interpolation of a soft resummation and a soft 
+ collinear one and constitutes a consistent description of a multi-scale process like semilep- 
tonic b ^ c decay; such interpolation is not unique and has been constructed on the basis of 
simplicity's requirements. The resummation formula in eq. (6) has been derived form "first 
principles" , namely from the universal properties of QCD radiation. Much attention has been 
paid to the correct inclusion of gluon-branching effects by means of a proper definition of gluon 
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transverse momentum. We have also presented in eq. (19) a simplified resummation formula 
which takes into account soft effects only and which can be applied as long as the charm quark 
is not too fast, as is the case in most of the phase space of semileptonic 6 — > c transitions. 



Resummed perturbation theory breaks down for final hadron masses as small as 




NP 



- ml 



'c ' 



(102) 



because of large soft effects controlled by the QCD coupling evaluated close to the Landau 
pole. Such effects can formally be factorized by the shape function for a massive quark in the 
final state [10]. In practice, one can model them with an effective ghost- less QCD couphng 
constructed on the basis of some analyticity requirement, along the lines of [26,27]. Our 
resummed calculation can be combined with the 0{as) triple differential distribution computed 
in [28, 29] by a proper matching procedure in order to obtain an uniform approximation in the 
whole phase space. 

We shortly mention here some phenomenological outputs of our results. Resummation ef- 
fects are expected to be largely sensitive to the charm mass and a phenomenological analysis 
could probably results in a more precise determination of m,.. Unlike the charmless channel, 
the b ^ c rate is known to be dominated by few hadronic states; given the high-quality data 
of the 5-factories, comparison with solid theoretical expectations could also provide a non- 
trivial check of local parton-hadron duality and, in general, it could give more information 
about hadron dynamics at the scale of few GeV's. The theoretical knowledge of 6 ^ c semilep- 
tonic spectra could also be of utihty for the inclusive extraction of the CKM matrix element 
\Vcb\- Furthermore in [26] a model was formulated which is in disagreement with the charm- 
less electron spectra below 2.2 GeV measured at the S-factories; a natural explanation is an 
under-subtracted charm background, which could be better estimated by means of our results. 

Our formula can also be used to describe beauty-mass effects in semileptonic top decays 
t b + W. A natural extension of our formalism would also allow a smooth description of 
heavy-flavor production in e'^e~ annihilation. 



as function of the center-of-mass energy Ec.m., from threshold, Ec.m. = '^mQ, up to higher 
energies, Ec.m. ^ ^mq. As far as soft effects are concerned this process may be seen as a 
crossed semileptonic decay. Coherence effects play a similar role as in the heavy flavor decay 
also in the production process (103): one produces a QQ pair in a color singlet state, which 
at the threshold, classically, sitting one on the top of the other, do screen their color and 
consequently emit no soft radiation. 



QQ, 



(103) 
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